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A new approach to the old problem of the predominance of prolate deformations among well 
deformed nuclei is proposed within the shell model framework. The parameter space is explored 
using the ensemble of random rotationally-invariant interactions. Subsets with rotational energy 
ratio E{A'^) / E{2'^) and the rigid-rotor relation between the quadrupole moment Q(2''") and the 
transition probability B{E2\2'^ — > 0^)) are found exhibiting prolate predominance. We identify 
matrix elements of the effective forces responsible for the predominance of prolate deformation. 



PACS numbers: 21.60.Cs, 21.10.Re, 21.60.Ka 

The majority of deformed nuclei have axially symmet- 
ric prolate deformation in their ground state (g.s.). The 
presence of rotational bands with typical energy intervals 
is the first signature of stable deformation. The states in 
a band are connected by strong quadrupole transitions 
obeying simple rigid rotor intensity rules. The same 
rules []| determine the expectation values of multipole 
operators. This is essentially a projection of nearly con- 
stant within the band intrinsic (body-fixed) quantities 
onto the laboratory (space-fixed) coordinate frame. The 
overwhelmingly positive sign of the intrinsic quadrupole 
moment manifests the dominance of prolate deformation. 
The Nilsson diagram shows a pattern of split spherical 
single-particle orbitals that seems to be more or less sym- 
metric with respect to the sign of deformation; the differ- 
ence in the liquid drop energy for the two signs of defor- 
mation is pretty small Q ■ While the results of the numer- 
ous mean field calculations for specific nuclei in general 
correctly predict the presence and the sign of deforma- 
tion, the underlying physical reason for the predominance 
of prolate shapes is unclear. 

It was argued by Lemmer |3(] that the kinetic energy 
should contain an additional contribution similar to the 
term —D1^ introduced by Nilsson f3| for interpolation 
between the harmonic oscillator and potential box. Be- 
ing split to time-conjugate pairs (m, — m) by deforma- 
tion, the large-€ spherical orbitals determine the single- 
particle (s.p.) occupancies in the way that makes the 
prolate case energetically favorable. Related arguments 
were given by Castel and Goeke 0], who showed that 
the collective inertial parameter is larger for the prolate 
deformation, and later by Castel, Rowe and Zamick Q 
with the help of self-consistency conditions There 
are also ideas based on the semiclassical analysis of the 
s.p. level density 0], simplest periodic orbits [§| and 
their bifurcations in a deformed cavity . The predomi- 



nance of prolate deformation analogously emerges in s.p. 
motion in metallic clusters [loj and for many biological 
objects, from molecular level to pollen grains [Hi. The 
latest analysis by Hamamoto and Mottelson [fj] starting 
with the statement that "the nature of the parameters 
responsible for the prolate dominance has not yet been 
adequately understood" proceeds along similar lines. De- 
tailed calculations and comparison between the harmonic 
oscillator potential and that of a spheroidal cavity show a 
different character of the mixing of spherical orbitals for 
the prolate and oblate cases. Being a surface effect, the 
difference should decrease in large systems. The authors 
of Ref. mention additional factors which were not 
sufficiently accounted for inclu ding the role of the spin- 
orbit potential, pairing effects [l^, and the presence of 
two kinds of particles. 

Below we suggest a new approach to this problem that 
was previously attacked from the viewpoint of the de- 
formed mean field and corresponding quasiparticle mo- 
tion. As the mean field itself is generated by the nucleon- 
nucleon interactions, it makes sense to take a step back 
and consider the many-body problem rather than the re- 
sulting s.p. motion. This would allow us to avoid ques- 
tions of self-consistency and take into account not only 
pairing and spin-orbit forces but the full inter-particle in- 
teraction. The natural framework is provided by the shell 
model (SM) with effective interactions; modern versions 
of the SM give good agreement with the data. The exact 
diagonalization of the Hamiltonian matrix, with all con- 
servation laws strictly respected, provides the eigenstates 
in the space-fixed frame; the spectrum and observables 
can be analyzed similarly to the experimental data. The 
problem requires exploration of the parameter space of 
the SM Hamiltonians in a sufficiently large orbital space. 

To explore the parameter space we will work with ran- 
dom interactions. All two-body matrix elements allowed 
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FIG. 1: The distribution of the Alaga ratio, Eq. for all cases selected by the energy ratio i?(4)/_E(2); the interaction 
strength is A = 0.5, panel a, and A = 4, panel b; the same for all cases with the sequence (0''',2"'"), panel c. 



by angular momentum, parity and isospin conservation 
are taken as random uncorrelated quantities. Among 
many realizations of such an ensemble we can identify 
those with definite rotational characteristics and study 
the abundance of such cases and their dynamic origin. 
Random interactions in many-body systems are usually 
analyzed for studies of quantum chaos iM, ,15i| in relation 
to the random matrix theory. In Ref. [16[ the ensemble 
was constructed with rotationally invariant interactions, 
and only the magnitudes of allowed matrix elements were 
taken randomly. The unexpected output was the discov- 
ery of the predominance, for an even number of particles, 
of the zero ground state spin, in spite of relatively low 
multiplicity of J = states in Hilbert space. This result, 
valid not only in a single- case but in the context of real- 
istic SM space 17] as well, was fully understood in terms 
of geometry of the parameter space only for the excep- 
tional case of the single j — 7/2 shell fis'l and for the sim- 
ilar results in the interacting boson model 19]. In more 
realistic fermionic models the full theory is still absent; 
the most promising are the ideas of geometric chaotic- 
ity of random angular momentum coupling l20l. l21l . The 
details were discussed in the review articles [22|, |23|, [24 ] . 

The rotationally invariant random interaction models 
provide a wealth of information. Any given set of random 
interaction parameters gives rise to a specific version of 
the SM. It was noticed, for example, that some sequences 
of the lowest excited states, like 0+ — 2+ — 4+ — appear 
with the enhanced probability. The energy ratios in such 
yrast sequences are spread around the values close to 
limiting cases of the collective models, such as vibrational 
and rotational bands [l^ 
with the lowest states 0+ 
the Alaga ratio, 
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In Ref. [22], the realizations 
- 2"*" were analyzed in terms of 



B{E2;2+ 0+)' 
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of the squared expectation value of the quadrupole mo- 
ment of the 2+ state to the reduced transition probability 
from this state to the g.s. The distribution of the Alaga 
ratio reveals two pronounced peaks, close to zero and 
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TABLE I: The probabilities of occurrence of collective prolate 
configurations (see text for details). 

at the value corresponding to the rigid rotor [l|]. This 
means that the stable mean fields generated by random 
interactions with high probability correspond either to 
near-spherical or to well deformed shape. The situation 
is similar in the case of the interacting boson model [l^ . 

We define the space-fixed quadrupole moment of the 
axially symmetric rotor with spin J [l| 



as 



Q{J) 



- j(j + 1) 

(J+ l)(2J + 3) 



J 



2j + r 



(2) 



where Qo is the intrinsic (body-fixed) quadrupole mo- 
ment, and we assume K — Q for the yrast band of an 
even-even nucleus. Prolate intrinsic shapes, Qo > 0, cor- 
respond to squeezed shapes around the axis of collective 
rotation, Q < 0. With the standard pj definition of the 
B(E2) probability for the yrast band of a well deformed 
rotor, the Alaga ratio ([!]) is equal to 4.10. 

In the single-j model [ll], where the spHtting of a 
j-level proceeds fan-like with approximate symmetry of 
oblate and prolate sides, the sign of Q will be random. 
Therefore we expect the crucial role of mixing of spheri- 
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cal orbitals. It was often suggested [20| that the spherical 
SM can reahstically describe deformation if one uses at 
least two s.p. orbits with Aj — 2. This approach suc- 
cessfully produces collective bands, such as in "^^Cr, and 
even in the double magic ^^Ni [27j. The collectivity was 
tested not only by the J( J+1) behavior, but also by large 
and consistent in band quadrupole moments and B{E2) 
strengths. As a generic example, we consider a model 
with four neutrons and four protons in the space of two 
spherical orbitals, O/7/2 and 1^3/2, of the same parity. 
This corresponds to the oversimplified SM description of 
"'^Cr, the nucleus with well known collective properties. 
Without interactions, we would have half-filled /7/2 or- 
bitals with zero quadrupole moment. In this space the 
most general rotationally invariant two-body interaction 
is described by 30 matrix elements. Table II below. We 
choose the random ensemble defined by the uniform dis- 
tribution of uncorrelated matrix elements between ~V 
and +V, and the only parameter is the ratio A = V/e, 
where e is the (non-random) spacing between the two or- 
bitals. With the statistics of 10000 realizations for each 
value of A, the results are presented in Table la. 

The case of very small A = 0.05 can be viewed as the 
single-j result. The column A^(0, 2) shows the number 
of realizations with the sequence 0"'',2+ of the g.s. and 
the first excited state. The probability of such sequences 
is strongly enhanced. Before applying additional require- 
ments we can find the sign of the expectation value ^ of 
the quadrupole moment Q(2+). As expected, the frac- 
tion N{Q < 0)/7V(0,2) of "prolate" cases, Q < 0, is 
stable near 50% (the third column) . Among these states 
there are cases with rotational properties. The constraint 
that the ratio E{'i+)/E{2+) be between 3 and 3.6 (the 
rigid rotor would give 3.33) leaves the number of states 
indicated in the column N{EA/ E2). In the column A'rot 
we count the states where, in addition to the energy crite- 
rion, the Alaga ratio ^ is between 3.90 and 4.30, an ar- 
bitrary but restrictive choice. The last column gives the 
fraction of states with negative sign of Q among all "ro- 
tational" states. The predominance of prolate configura- 
tions is practically independent on the relative strength 
of mixing between the orbitals if this strength is sufficient 
for the onset of deformation. The full histogram for the 
distribution of the Alaga ratio is shown in Fig. 1 for 
two values of the interaction strength; this figure shows 
also the statistics without any rotational selection cuts. 
The prolate deformed peak is quite narrow. With the 
only energy criterion (0+, 2+) applied, we see also a large 
peak for spherical configurations with A close to zero, as 
well as intermediate background situations. In this large 
set, the prolate cases appear only in slightly more than 
50% cases. We can also note that at very weak mixing, 
A = 0.05, the numbers iV(0, 2) and especially N{EA/E2) 
fall sharply, and the number of rotational configurations 
satisfying our criteria is very small. This emphasizes the 
importance of orbital mixing. 
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FIG. 2: The average matrix elements of the two-body inter- 
action. The matrix element number is as in Table II. 

Going from the isospin-symmetric case with N = Z 
to asymmetric one, we consider next the system with six 
neutrons and four protons assuming the same set of two 
s.p. levels. The results are shown by Table lb. The frac- 
tion of states with Q < among all sequences (0"'",2"'") 
is still close to 1/2. The total number of rotational cases 
is now significantly lower. At the same time, the frac- 
tion of prolate configurations is even higher than in the 
symmetric case. The results are essentially the same in 
the N = Z case for the inverted level scheme when the 
spherical P3/2 level is placed below the /7/2 level, see Ta- 
ble Ic. Again for the very weak mixing, A = 0.01, we 
do not find any rotational states; in this case all eight 
nucleons occupy the lowest P3/2 level creating the closed 
shell. With strong level mixing, 45% of all states be- 
long to the (0"'",2"'") sequence. Clearly, the occurrence 
of prolate deformation is even higher than for the nor- 
mal fp level sequence. This confirms the presence of 
the effect we mentioned in the introduction, namely the 
influence of the splitting of the level with higher £ when 
1ow-|to| components steeply go down with the prolate de- 
formation. The last case we consider here contains four 
protons and four neutrons in the space of two levels of 
opposite parity, f^^2 and 59/2, Table Id. Here the fraction 
of the sequences 0^,2+ is lower, and practically all ro- 
tational states have oblate deformation. It is clear that, 
due to parity conservation, direct mixing of split levels by 
the deformed mean field is impossible and only pairwise 
transfers of particles can play a role. 

The main new dimension brought in by our approach 
is in the visualization of the entire parameter space. The 
mechanism of the prolate predominance is revealed as we 
identify the realizations of the ensemble responsible for 
the effect. Fig. 2 shows all 30 reduced matrix elements 
{jij2\\V\\j3j4){JT) allowed by angular momentum and 
isospin for the strong coupling limit, A = 4.0, in case (a) 
of Table II. The average value of each matrix element over 
the whole ensemble is zero, with small fluctuations (solid 
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TABLE II: List of average matrix elements for the case of 
Table la and A = 4.0. 

line). The dashed line shows the average over prolate 
rotational samples. Table II displays the ordering of the 
matrix elements and their average numerical values. Al- 
though all matrix elements are enhanced compared to the 
level of fluctuations, one can conclude that certain ma- 
trix elements are crucial for the transition from spherical 
shape to predominantly prolate axial deformation. We 
see the exceptional role of amplitudes 9-12 describing the 
transfer of a single nucleon fi^2 ^ Vzji the interac- 
tion with another /7/2 nucleon, regardless of its isospin. 
Such a process is forbidden for the orbitals of opposite 
parity. This agrees with the idea discussed in Ref. 
The monopole pairing given by amplitudes 5, J 5 and 29 is 
not effective in this transition. Contrary to that, we see 
a large negative amplitude of quadrupole pair transfer 
IQ. Large quadrupole-quadrupole forces in the particle- 
hole channel correspond to matrix elements 20-24 in the 
particle-particle channel. They induce the collectiviza- 
tion process and formation of the proper symmetry of 
the mean field after mixing p- and /-orbitals. 

In summary, for the first time we performed the explo- 
ration of the parameter space that serves as an arena for 
effective nucleon-nucleon interactions building the sta- 
ble deformed mean field. We show the decisive role of 
the mixing between the valence spherical orbitals of the 
same parity split as a function of deformation. This mix- 



ing, different for the two sides of the axial deformation 
diagram, makes the prolate deformation energetically fa- 
vorable with high probability. This picture is supported 
by the statistical analysis of the random interaction en- 
semble and by singling out the responsible components 
of the effective interaction. The process is amplified by 
the presence of two kinds of nucleons. Although we have 
shown in detail in this Letter only cases of simple config- 
urations, the results are quite generic for a small system. 
The effect is driven by the surface orbitals but it still re- 
mains to see if it indeed d isap pears in the macroscopic 
limit as predicted in Ref. [l^. The random interaction 
ensemble provides a new powerful tool for understanding 
the many-body mechanisms of collective phenomena. 
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